arXiV: astro-ph/0112406v 1 17 Dec 2001 


RAPID FORMATION OF ICE GIANT PLANETS 


Alan P. Boss, George W. Wetherill, and Nader Haghighiponr 


Department of Terrestrial Magnetism 
Carnegie Institution of Washington 
5241 Broad Branch Road, N.W., Washington, D.C. 20015-1305 
(202)478-8858, FAX (202)478-8821, boss@dtni.ciw.edu 


Number of manuscript pages: 20 
Number of figures: 0 
Number of tables: 1 

[Submitted as a Note to Icarus] 


1 


Suggested running head: Ice Giant Planet Formation 


Send correspondence to: 


Dr. Alan P. Boss 

Department of Terrestrial Magnetism 
Carnegie Institution of Washington 
5241 Broad Branch Road, N.W. 
Washington, D.C. 20015-1305 


2 



ABSTRACT 


The existence of Uranus and Neptune presents severe difficulties for the core accretion 
model for the formation of ice giant planets. We suggest an alternative mechanism, namely 
disk instability leading to the formation of gas giant protoplanets, coagulation and settling 
of dust grains to form ice/rock cores at their centers, and photoevaporation of their gaseous 
envelopes by a nearby OB star, as a possible means of forming ice giant planets. 
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INTRODUCTION 


The robustness of the current theory of the origin of the Solar System deteriorates 
sharply with increasing distance from the Sun. The innermost, terrestrial planets are 
widely believed to have formed by the collisional accumulation of rocky planetesimals, 
culminating after ~ 50 Myr in giant impacts between planetary-sized bodies (Canup and 
Righter 2000). The gas giant planets, Jupiter and Saturn, were formed by core accretion 
in the conventional view, where a ~ 10 Earth-mass (~ IOM 0 ) solid core forms hrst by 
collisional accumulation, and then accretes a gaseous envelope from the solar nebula, a 
process requiring an estimated 8 Myr in the standard model (Pollack et al. 1996). However, 
most young stars lose their gas disks in a few Myr or less in regions of low-mass star 
formation (Briceho et al. 2001) and on even shorter time scales in OB associations (Bally et 
al. 1998), meaning that core accretion can only produce gas giant planets in the rare long- 
lived disk. Alternatively, gas giant protoplanets might form rapidly (in ~ 10^ yrs) by a disk 
instability, where the gas in a marginally gravitationally unstable disk forms clumps which 
then contract to planetary densities (Boss 1998, 2000, 2001a). Disk instability operates 
fast enough to form gas giant protoplanets in even the shortest-lived protoplanetary disks. 

The problems facing the formation in situ of the ice giant planets, Uranus and Nep¬ 
tune, by the core accretion mechanism are even more severe than those facing the formation 
of Jupiter and Saturn. If the ice giants formed in situ, the lower surface density of solids 
and longer orbital periods require that collisional accumulation proceed much slower than 
in the gas giant planet region (inside ~ 10 AU). Even more importantly, because the escape 
velocity from the Solar System at 20 AU to 30 AU (the approximate orbital distances of 
Uranus and Neptune, respectively) is ~ 8 km s“^, comparable to orbital velocities and 
to the relative velocities between growing embryos, the effect of mutual encounters is to 
excite orbital eccentricities so much that the embryos cross the orbits of Saturn or Jupiter, 
and thus are either ejected on hyperbolic orbits, lost by impact with the gas giant planets, 
or kicked into cometary orbits (Lissauer et al. 1995). Ice giant planets cannot form in 
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the standard model (Levison and Stewart 2001). Possible means for salvaging the colli- 
sional accretion mechanism inclnde (Levison and Stewart 2001) invoking some sort of drag 
force to damp orbital eccentricities, or rnnaway accretion of a single embryo to Uranns-size 
(Bryden et al. 2000) rather than oligarchic growth of mnltiple embryos (Koknbo and Ida 
1998). Perhaps the best snggestion for forming the ice giant planets by collisional accretion 
is the proposal that the ice giants were formed between Jnpiter and Satnrn and then were 
gravitationally scattered ontward (Thommes et al. 1999). However, the snccess of this 
solntion depends largely on the ability of core accretion to form gas giant planets prior to 
the removal of the disk gas, which seems to be limited, considering what we know abont 
gas disk lifetimes. 

Here we propose a completely different mechanism for the formation of the ice giant 
planets, motivated by recent theoretical work on gas giant protoplanet formation by the 
disk instability mechanism (Boss 1998, 2000, 2001a) and by observations of the erosion of 
Orion Nebnla protoplanetary disks by extreme nltraviolet (EUV) irradiation (Bally et al. 
1998). In a nntshell, we propose forming the ice giant planets rapidly throngh the sedi¬ 
mentation of solids in massive gas giant protoplanets formed by disk instability, followed 
by removal of most of their gaseons envelopes by intense EUV irradiation from a nearby 
massive star. We believe that this mechanism may represent a snperior means for form¬ 
ing ice (and gas) giant planets in the majority of protoplanetary disks. Terrestrial planet 
formation is still likely to occnr even when gas giant planets are formed qnickly by a disk 
instability (Kortenkamp and Wetherill 2000), and may even be facilitated (Kortenkamp et 
al. 2001). 


DISK INSTABILITY 

The disk instability mechanism is capable (Boss 2000) of forming mnltiple-Jnpiter- 
mass clnmps in a marginally gravitationally nnstable disk with a snrface density at 5 AU 
comparable to that believed necessary to form Jnpiter by core accretion. Extremely high 
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spatial resolution is necessary in order for self-gravitating clumps to form and to persist; 
three dimensional (3D) calculations with at least 10® grid points seem to be necessary 
(Boss 2000), which should be compared to the 25,000 particles used by Laughlin and 
Bodenheimer (1994) in their 3D smoothed particle hydrodynamics (SPH) disk simulation. 
Disk instability appears to be capable of forming self-gravitating clumps even when a 
rigorous thermodynamical description (energy equation, radiative transfer in the diffusion 
approximation, detailed equations of state, and dust grain opacities) of the disk is included 
(Boss 2001a), because the disk cooling time is comparable to the relevant dynamical time 
scale, namely the orbital period. Models where artihcial viscosity is used to heat the gas, 
on the other hand, tend to damp the growth of nonaxisymmetry (Nelson et al. 2000, 
Pickett et al. 2000a,b). Artihcial viscosity was not employed in the 3D models of Boss 
(1998, 2000, 2001a). 

In spite of these promising results, it must be admitted that the present disk insta¬ 
bility models cannot yet be considered dehnitive, because of inherent limitations, such 
as the spatial resolution employed even in the highest spatial resolution models to date 
(Boss 2000, 2001a), and doubts regarding the long-term survival of the clumps. However, 
models currently in progress by one of us (APB) suggest that with specialized techniques 
to improve the spatial resolution of, e.g., the Poisson solver for the gravitational potential, 
the disk instability mechanism becomes increasingly robust. 

While the Boss (2000, 2001a) models were limited to clump formation inside 20 AU, 
disk instability models with greater radial extent imply that clumps could still form at 
somewhat greater distances in a suitable disk (Nelson et al. 1998; see Boss 2001a for a 
discussion of the validity of the “locally isothermal” approximation used by Nelson et al. 
1998 and Boss 2000, compared to the approximations used in modeling heating and cooling 
processes by Nelson et al. 2000). The type of protoplanetary disk that would be necessary 
to produce clumps at 20 AU to 30 AU, i.e., a disk with a gravitational stability parameter 
Q < 1.5 in that region, would have a total mass inside 30 AU of ~ 0.13 solar masses (Mq), 
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with ~ 0.04Mq residing between 20 AU and 30 AU, based on extending previons disk 
instability models (Boss 2000, 2001a) to a radins of 30 AU. The assnmed resnlt of snch 
a disk model, fonr gas giant protoplanets spaced from ~ 5 AU to ~ 30 AU, is plansible 
based on the previons modeling, bnt has not as yet been calcnlated convincingly. 

CORE FORMATION 

Table 1 lists the cnrrent estimates (Gnillot 1999, Hnbbard 1992) for the total, core, 
and envelope masses for the gas and ice giant planets, where the gaseons envelope mass is 
assnmed to be the difference between the total mass and the inferred core mass. Note that 
a range of core masses is possible for Jnpiter and Satnrn, and the midpoint of this range is 
nsed in Table 1; the core masses of Uranns and Neptnne are similarly nncertain, so nominal 
valnes are employed. Also shown in Table 1 is the mass of the gas giant protoplanet that 
wonld be needed to acconnt for the inferred core masses, assnming a gas to dnst mass ratio 
of 50:1 and that all the dnst grains coagnlated together and sedimented to the center of the 
protoplanet. The hnal row in Table 1 then lists the amonnt of hydrogen and helinm gas 
that wonld need to be removed to leave behind a planet with the observed total mass. We 
do not consider here the qnestion of the heavy metal enrichment of the envelopes of the gas 
giant planets, as for either core accretion or disk instability, this envelope enrichment is 
likely to be associated at least in part with the ingestion of planetesimals after the planet 
has formed. 

Insert Table 1 

Clnmp formation can occnr very rapidly in a disk instability, on the dynamical time 
scale of the orbital period, or within abont a few hnndred years at an orbital distance of 
25 AU. Once a well-dehned clnmp forms, the dnst grains within the clnmp will begin to 
sediment toward the center of the clnmp. Dnst grain growth and sedimentation to the 
disk midplane may have already begnn at the time when the clnmps formed, giving a head 
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start to this process, but to be conservative, we assume that no growth has yet occurred 
when the clumps form. Using the standard approach for calculating the growth of dust 
grains by differential settling in the solar nebula (Weidenschilling 1988), specialized to the 
case of a spherical protoplanet, we obtain a characteristic growth time (Boss 1998) of 

_ _ Cs 

“ TifppGR^ 

where Cg — Cs(T) is the sound speed, / is the dust to gas mass ratio, Pp is the density of 
the protoplanet, G is the gravitational constant, and R is the radius of the protoplanet. 
For nominal values at 25 AU of the gas temperature T = 50iF, / = 0.02, pp = 10“® g 
cm“^, and R = 2 AU, we obtain Ta ~ 10 yrs. The time to grow from size uq to a is 
~ 3Taln{a/ao), so growing from ag = 0.1 pm to 1 cm-size requires ~ 10^ yrs. Cm-sized 
particles can then sediment to the center of the protoplanet in a time 

3cs 

47raGpa ’ 

where pa ~ 3 g cm“^ is the density of a dust grain. For the nominal values, the sedimen¬ 
tation time of cm-sized grains is Tc ~ 2 x 10^ yrs. Hence the solids inside a gas protoplanet 
at ~ 25 AU should coagulate and form a central core on a time scale of ~ 3 x 10^ yrs. 

This estimate of the time scale for the formation of a solid core assumes that the 
protoplanet is in radiative equilibrium rather than in convective equilibrium. Wuchterl 
et al. (2000) calculated that a gas giant protoplanet, starting from conditions which 
approximately characterize its formation during a disk instability, is likely to become nearly 
fully convective within about 100 yrs. Turbulent motions driven by convection may thus 
need to be taken into account in assessing whether or not dust grains can settle to the 
center of a protoplanet. While this situation is clearly deserving of further study, we can 
suggest two possible outcomes here. First, considerable dust grain coagulation in the solar 
nebula may have already occurred by the time at which the disk instability occurs, so 
the assumption of starting from 0.1 pm-size may be too severe. All that is required is 



for the solids to be small enough to remain coupled to the gas sufhciently that they are 
incorporated into the gaseous clumps during the disk instability phase. If, for example, 
10-cm-sized solids had formed prior to the instability, they would settle to the center of 
the protoplanet in a few 100 yrs, perhaps prior to the onset of widespread convection. 
Alternatively, detailed calculations of the coagulation and vertical settling of solids in the 
solar nebula show that even in the presence of moderate turbulence (e.g., a turbulent 
velocity of 0.01 km s“^), settling to the midplane at 1 AU can occur in about 2000 yrs 
(see Fig. 9 of Weidenschilling and Cuzzi 1993), i.e., in about the same time as for a non- 
turbulent nebula (see their Fig. 6). These calculations suggest that even a fully convective 
protoplanet need not present an insurmountable barrier to core formation. Afterall, for 
every updraft there is a downdraft, in isotropic turbulence. It remains to be seen if the core 
formation process can produce a core which is nearly hydrogen-free, as would be required 
to produce Uranus and Neptune, and if a hydrogen- and helium-rich envelope will be stable 
during and after this process. Detailed radiative hydro dynamical models will be necessary 
to address these points (e.g., Wuchterl 1993). 

Water ice is stable for temperatures less than 180 K (Pollack et al. 1994) at a density of 
10“® g cm“^. By comparison, a newly-formed gas giant protoplanet may have a maximum 
temperature of about 130 K at a density of 10“® g cm“^ (see Figs. 2c,d of Boss 2001a), 
suggesting that water will be in a solid form, at least initially. The question of whether 
water and other ices will remain solid during the core formation phase remains to be 
investigated in detail. If the ices sublimate, icy core formation could be halted (the critical 
temperature for water ice is 647 K at a pressure of 218 atm). Addressing this issue will 
require a detailed two-fluid calculation of the heating and cooling processes in a protoplanet 
during the core formation process. At any rate, more refractory species should continue 
to settle to the core; Uranus’s core is thought to be composed of 2/3 ice and 1/3 rock 
(Hubbard 1997). 
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PHOTOEVAPORATION 


An isolated giant gaseons protoplanet shonld contract toward planetary densities on 
a time scale of ~ 10^ yrs (Bodenheimer et al. 1980, Cameron et al. 1982). We hypothesize 
that dnring this period, a massive star formed in the same star-forming region as the solar 
nebnla, and began to prodnce copions EUV irradiation, similar to the sitnation in the Orion 
Nebnla clnster today. OB associations like Orion appear to be the birthplace of most stars 
(Bally et al. 1998), so an Orion-like environment shonld have inflnenced the formation 
of many planetary systems, possibly inclnding onr own. Signihcant losses (Richling and 
Yorke 1998) of disk gas can occnr in as little as ~ 10^ yrs for the protoplanetary disks 
(proplyds) close to the Orion Trapezinm star 0^ C Ori, while disk snrvival times are 
estimated (Bally et al. 1998) to be of order 10^ yrs, depending on the disk mass. For a 
disk massive enongh to nndergo a gravitational instability, the snrvival time may be ~ 10® 
yrs (see below). Note that we have assnmed here that the disk instability which forms 
the giant protoplanets occnrs prior to the onset of EUV irradiation, bnt this need not 
necessarily be the case. Becanse external EUV radiation photoevaporates a disk from the 
ontside in, it is conceivable that at the optically thick disk midplane, the disk may remain 
snfficiently cool (25 K to 50 K) for a disk instability to proceed even after irradiation has 
commenced. 

EUV irradiation, however, cannot evaporate disk gas closer than abont 5 to 10 AU 
from a solar-mass protostar, becanse the gas temperatnres prodnced ( 10^ K) are insnfh- 

cient to permit a thermally-driven wind to exit the gravitational potential well of the star 
(p. 761 of Johnstone et al. 1998). This hot gas will remain gravitationally bonnd to the 
star, forming a protective, ronghly spherical halo with a radins of abont 5 to 10 AU. Hence 
we expect that protoJnpiter will remain shielded from the EUV flnx by the persistence of 
the inner solar nebnla inside this halo, protoSatnrn will be only partially shielded by the 
halo, while protoUranns and protoNeptnne, orbiting ontside the halo, will soon become 
nnshielded as the disk evaporates from the ontside in, and so they will begin to lose their 
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gaseous envelopes. The proto-ice giants evidently need to lose ~ 2 Jupiter masses (Mj) of 
gas (Table 1), with the amount needed to be removed decreasing with decreasing orbital 
distance and falling to zero inside 10 AU. The mass loss rate for gas removal from disks or 
clumps close (~ 10^^ cm) to 0^ C Ori by EUV irradiation is estimated (Johnstone et al. 
1998) to be 


M ~ 8 X 10 Tc Mqut 

where Tc is the radius of the clump in cm. The radii of giant gaseous protoplanets are 
limited by tidal truncation by the protosun, with the critical tidal radius being propor¬ 
tional to orbital distance (Boss 1998), so disk instability models (Boss 2001a) showing the 
formation of clumps with radii of ~ 0.6 AU at 10 AU would predict clumps with radii of ~ 
1.5 AU at a 25 AU distance. Such a clump would then lose mass at a rate of ~ 2 x 10 ~^Mj 
yr“^, sufficient to lose the required amount of envelope gas in ~ 1 Myr. This mass loss 
rate may also be used to estimate the time scale for EUV radiation to hrst remove the 
disk gas (~ 0.04Mq) residing between 20 AU and 30 AU in a marginally gravitationally 
unstable disk, which is ~ 1 Myr. 

While an isolated clump may contract to considerably smaller radii than 1.5 AU in 
about 10® yrs (Bodenheimer et al. 1980, Cameron et al. 1982), once EUV radiation reaches 
the protoplanet, it will be embedded in a high temperature (10,000 K) plasma that may 
slow the escape of its radiation, stop any further contraction, and may even expand the 
protoplanet as its outer layers are being stripped away (Cameron et al. 1982). In fact, once 
gas giant protoplanets form, they are likely to accrete or otherwise clear the disk gas in their 
vicinity, opening up disk gaps, though this phase of evolution has not been encountered 
in the simulations to date. Unless the OB star lies exactly in the disk plane and EUV 
light cannot be scattered onto the protoplanets, the star’s EUV radiation should begin to 
strip the outer layers of the protoplanets immediately, well before the unaccreted disk gas 
has been photoevaporated. Hence the estimated time scale of ~ 1 Myr for protoplanet 
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envelope removal may still be approximately correct, though it remains to be seen whether 
the onset of photoevaporation in the outermost layers can signihcantly slow the contraction 
of a protoplanet. If not, then the OB star might have to be located closer to the disk than 
is assumed here, nominally 10^^ cm, in order to have envelope removal occur on a time 
scale comparable to that for protoplanet contraction. 

On the other hand, it must be noted that the protoplanet contraction time estimates 
of Bodenheimer et al. (1980) and Cameron et al. (1982) are based on spherically symmet¬ 
ric models, where rotation is ignored. In reality, gas giant protoplanets will be rotating 
very rapidly as they form (Boss 2001a), and their rate of contraction will be limited not 
only by the rate at which they can radiate the compressional energy produced by con¬ 
traction, but also by the rate at which they can transfer angular momentum out of their 
central regions. These protoplanets should thus look like rotationally-flattened polytropes, 
with outer regions where centrifugal forces provide signihcant support against self-gravity. 
Because when the protoplanets hrst form, they are only marginally self-gravitating (Boss 
2001a), and because of the heating associated with further contraction toward planetary 
densities, these protoplanets should have gravitational instability parameters Q > 1.5 dur¬ 
ing their evolution, and thus should not experience sub-fragmentation. It is unclear what 
processes will control the outward transport of angular momentum in these objects, but 
it is certain that the time scale for contraction must be longer than that in the absence of 
rotation. Hence contraction times based on spherically symmetric models should be lower 
limits, and the true contraction times should permit signihcantly more time for EUV to 
photoevaporate the rotationally distended envelopes of the outer protoplanets. Detailed 
modeling of the evolution of a rotating gaseous protoplanet is needed to address this key 
issue. 

EUV radiation will not be able to remove the entire gaseous envelope of protoUranus 
and protoNeptune, for the same reason that the solar nebula inside 5 to 10 AU cannot be 
removed. The inner radius Vg of EUV removal (Johnstone et al. 1998) is 
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where Mp is the planet mass and Cg = 10 km s“^ for gas heated by EUV radiation. For a 
hnal planet mass of ~ Vg — 3 x 10^ cm, slightly larger than the present-day radii 

of Uranns and Neptnne of ~ 2.5 x 10® cm, implying that EUV envelope stripping might 
stop short of the solid core, leaving behind a small gaseons envelope, as is observed. Other 
choices of Mp wonld yield other valnes of Tg] this analysis simply shows that Mp ~ 15M0 
is consistent with EUV irradiation and the present radii of the ice giant planets. Satnrn 
appears to be an intermediate case between the envelope losses snffered by the ice giants 
and the envelope retention enjoyed by Jnpiter, dne to its location close the point where 
the solar nebnla conld no longer be evaporated by external EUV radiation. 

Removal of most of the gaseons envelopes of protoUranns and protoNeptnne thns 
reqnires ~ 1 Myr of snstained EUV irradiation from a nearby massive star. Is this likely? 
The lifetimes of massive (> 8 Mq) stars can exceed ~ 10 Myr, and the ionizing stars of 
the Trapezinm clnster (Prosser et al. 1994) have already lived for 0.3 to 2 Myr. Star¬ 
forming OB associations (Elmegreen et al. 2000) have lifetimes on the order of 10 Myr, 
so ~ 1 Myr of snstained EUV irradiation from an external sonrce seems qnite possible. 
The solar nebnla need not have remained at a hxed distance from a star similar to 0^ 
C Ori dnring EUV irradiation: Scally and Clarke (2001) have shown that stars orbiting 
thronghont the Orion Nebnla clnster will lose considerable disk mass (0.01 to 1 Mq) by 
UV photoevaporation within abont 2 Myr. The early Snn itself (Shn et al. 1993) may 
have prodnced enongh EUV radiation to remove the onter solar nebnla gas beyond ~ 9 
AU over a time of ~ 10^ yrs, bnt protoUranns and protoNeptnne wonld have contracted to 
planetary-size dnring this period, making EUV evaporation inefficient. A strong external 
sonrce of EUV radiation seems to be reqnired for onr scenario to work. 
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CONCLUSIONS 


While we have focused on EUV ionizing radiation (A < 912A), far-ultraviolet (FUV) 
radiation (A > 912A) can also play an important role in evaporating protoplanetary disks 
(Johnstone et al. 1998). The four luminous Trapezium stars provide the EUV and FUV 
photons which bathe the protostars in the Orion Nebula cluster. FUV radiation will aid 
in the removal of gaseous envelopes from the proto-ice giants. 

Because most stars form in OB associations, EUV and FUV irradiation may play a 
much more important role in planetary system formation than has previously been thought 
to be the case. As we have shown here, EUV and FUV irradiation need not prevent the 
formation of gas or ice giant planets. Given that the emerging results (Butler et al. 2000, 
Boss 2001b) from the ongoing census of extrasolar gas giant planets give no indications that 
gas giant planets are rare, it seems that star-forming regions like Orion must be capable 
of forming at least gas giant planets, and presumably entire Solar Systems. 

APPENDIX: ORBITAL STABILITY PRIOR TO IRRADIATION 

Considering that the giant protoplanets required in Table 1 are more massive than the 
existing giant planets, one question which arises regards the stability of their orbits - can 
these protoplanets remain on stable orbits long enough for EUV irradiation to reduce their 
masses to their current values? The orbital stability analysis of Chambers et al. (1996) sug¬ 
gests that systems of more than two planets (with masses mi, m2) will be orbitally unstable 
around a 1 Mq star if their semimajor axes (ai,a 2 ) are separated by less than about 10 
Rh, where the mutual Hill radius is dehned to be Rh = [(wi -|-m2)/(3MQ)]^/^(ai -|-a2)/2. 
For the four giant planets shown in Table 1, located at the current orbital distances of 
the outer planets, all three pairs of planets are located closer together than 10 Rh (5.9 
Rh, 6.1 Rh, and 3.9 Rh, respectively, moving outward), implying instability. However, 
the time scale for the system to disrupt is estimated to be of the order of 10 Myr for the 
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innermost three giant planets, with separations of ~ QRr (extrapolating the resnlts shown 
in Fignre 3 of Chambers et al. 1996 to three planets, each with a mass 10 ^Mq), thongh 
the ontermost planet has a smaller separation (3.9 Rr) from its adjacent planet, which 
shonld shorten this time scale. In order to fnrther test the stability, one of ns (NH) has 
nsed the SWIFT symplectic integrator code of Levison and Dnncan (1994) to integrate the 
three dimensional orbits of fonr planets in orbit aronnd the Snn, starting with the cnrrent 
positions of the fonr onter planets (CHO: Cohen et al. 1973) and with the masses given in 
Table 1. With the CHO initial conditions, the system is stable for at least 10 Myr; only rel¬ 
atively small changes (~ 10%) in semimajor axes are evident dnring this period. For other 
arbitrary choices of the initial conditions, however, the time scale for orbital instability 
can be as short as a few 0.1 Myr. The SWIFT calcnlations do not inclnde the interactions 
of the planets with the disk gas, and so cannot be considered dehnitive - for example, the 
protoplanets may have formed with larger mntnal separations and then nndergone orbital 
migration dne to disk interactions, prior to or dnring photoevaporation of the onter disk 
gas. While the orbital stability of the nominal initial conhgnration implied by Table 1 is 
an important qnestion to consider in fntnre work, we believe that the SWIFT integration 
from the CHO initial conditions snggests that orbital stability for several Myr or more is 
at least plansible, which is a time period snfhciently long for EUV irradiation to lower the 
ice giant protoplanet masses to their present valnes. 
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TABLE 1 


Total, inferred core (Guillot 1999, Hubbard 1992) and resulting envelope masses (in Earth 
masses, M^) for the gas and ice giant planets, masses of the gas giant protoplanets (GGPP, 
in Jupiter masses, Mj) needed to form the observed cores by the disk instability mecha¬ 
nism, and masses of envelope gas that must be removed by EUV irradiation (in Mj) to 
yield the observed planet. 


Mass 

Jupiter 

Saturn 

Uranus 

Neptune 

Total (M 0 ) 

318 

95 

15 

17 

Gore (M 0 ) 

6 

12 

12 

16 

Envelope {M^) 

312 

83 

3 

1 

GGPP (Mj) 

1 

2 

2 

2.5 

Removed (Mj) 

0 

1.7 

1.95 

2.45 
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